Hints on halo evolution in SFDM models with galaxy observations 
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A massive, self-interacting scalar field has been considered as a possible candidate for the dark 
matter in the universe. We present an observational constraint to the model arising from strong 
lensing observations in galaxies. The result points to a discrepancy in the properties of scalar field 
dark matter halos for dwarf and lens galaxies, mainly because halo parameters are directly related 
to physical quantities in the model. This is an important indication that it becomes necessary to 
have a better understanding of halo evolution in scalar field dark matter models, where the presence 
of baryons can play an important role. 

PACS numbers: 95.30.Sf, 95.35.+d, 98.62.Gq, 98.62.Sb. 



I. INTRODUCTION 

The nature of dark matter (DM) remains elusive today, 
even though a generic cold particle weakly coupled to the 
standard model seems to be the most promising candi- 
date [T] . Treating DM as a bunch of classical particles is 
an appropriate effective description for many physical sit- 
uations. However, if DM is composed of bosons, the zero 
mode can develop a non- vanishing expectation value; this 
effect is usually known as Bose-Einstein condensation. A 
condensed phase does not admit a description in terms 
of classical particles, and the concept of a coherent exci- 
tation (i.e. a classical field) is more appropriate for prac- 
tical purposes [2J. A specific realization of this scenario 
can be provided by the axion [3] , see also [I] . 

In this paper we shall explore the lensing properties 
of a generic model of DM particles in a condensate, and 
compare the conditions necessary to produce strong lens- 
ing with those required to explain the dynamics of dwarf 
galaxies. As a result we will get some insight into halo 
evolution arising from this type of models. 

In particular, we will consider the case of a com- 
plex, massive, self-interacting scalar field <j) satisfying 
the Klein-Gordon (KG) equation, Dd> — (mc/K) 2 cb — 
\\<t>\ 2 <t> = 0, with the box denoting the d'Alembertian 
operator in four dimensions. For those natural situa- 
tions in which the scalar field mass m is much smaller 
than the Planck scale, mpi anc k = (hc/G) 1/2 , such that 
A = Amp lanck /47rm 2 » 1, the coherent (self-gravitating, 
spherically symmetric) solutions to the KG equation ad- 
mit a very simple expression for the mass density [6] , 
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As usual we will refer to this m odel as scalar field dark 
matter (SFDM). Here r max = y / n 2 A/2 (h/mc) is a con- 
stant with dimensions of length (notice that r max is just 



the Compton wavelength of the scalar particle, h/mc, 
scaled by a factor of order A 1 / 2 ), and p c the density at 
the center of the configuration. The mass density profile 
in Eq. ([I]) leads to compact objects of size r max , and typ- 
ical masses, 4p c r 3 nax /7r, that vary from configuration to 
configuration according to the value of the central den- 
sity. 

Eq. was obtained without taking into account the 
gravitational influence of any other matter sources, and 
assuming that all the scalar particles are in the conden- 
sate. It has been used as a first order approximation to 
describe the distribution of matter in dwarf shperoidals, 
which are expected to be DM dominated [6 8J . The mass 
distribution would be smooth close to the center of these 
galaxies, alleviating the cusp/core problem motivated by 
the discrepancies between the observed high resolution 
rotation curves and the profiles suggested by N-body sim- 
ulations [S]; see however [TU] , 

The dynamics of dwarf galaxies suggests a self- 
interacting scalar with m 4 /A ~ 50 — 75(eV/c 2 ) 4 , (i.e. 
''max ~ 5.5 — 7Kpc), and typical central densities of the 
order of p c ~ lO _3 M0/pc 3 , see Ref. We are aware 
that Milky Way size galaxies are, at least, an order of 
magnitude larger than this value of r max , and then they 
do not fit in this model as it stands. Nonetheless, if 
not all the DM particles are in the condensate, there is 
a possibility to have gravitational configurations where 
the inner regions are still described by the mass density 
profile in Eq. ([!]) , wrapped in a cloud of non-condensed 
particles |llj . For the purpose of this paper we do not 
need to specify the complete halo model. This is because 
strong lensing is not very sensitive to the mass distribu- 
tion outside the Einstein radius, at most of the order of 
a few Kpc, just bellow the expected value of r max . Wc 
could not neglect the exterior profile of the halo if we were 
interested, for instance, in weak lensing observations. 
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II. LENSING PROPERTIES OF SFDM HALOS 

In the weak field limit the gravitational lensing pro- 
duced by a mass distribution can be read directly from 
the density profile. As usual we assume spherical symme- 
try, and use the thin lens approximation, that is, the size 
of the object is negligible when compared to the other 
length scales in the configuration, i.e. the (angular) dis- 
tances between the observer and the lens, -Dol, the lens 
and the source, -DlSj and from the observer to the source, 
Dos- 

Under these assumptions the lens equation takes the 
form 
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with p and denoting the actual (unobservable) angular 
position of the source, and the apparent (observable) an- 
gular position of the image, respectively, both measured 
with respect to the line-of-sight [12]. The (projected) 
mass enclosed in a circle of radius £, M(£), is defined 
from the (projected) surface mass density, £(£), through 

dzp(z,$, M(£) = 2tt / deem')- (3) 

-oo JO 

Here £ = Dq^Q is a radial coordinate in the lens plane, 
and z a coordinate in the orthogonal direction. Finally 
S cr = c 2 Dos / '47rG DolDls is a critical value for the sur- 
face density. 

In general, Eq. ^ will be non-linear in 0, and it could 
be possible that for a given position of the source, /3, there 
would be multiple solutions (i.e. multiple images) for 
the angle 0. This is what happens in the strong lensing 
regime to be discussed below. One particular case is that 
with a perfect alignment between the source and the lens, 
that actually defines the Einstein ring, with an angular 
radius of E = 6(fi = 0). 

For a SFDM halo, and in terms of the normalized 
lengths e = £/r max and = z/r max , the surface mass 
density takes the form 



^SFDM 
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dz* , (4) 



with < e < 1 and z max = \A — £* • A similar 
expression can be obtained for the mass enclosed in a 
circle or radius £, see Eq. Q above. Here we are not 
considering the effect of a scalar cloud surrounding the 
condensate. For r < r max this will appear as a projection 
effect, which is usually considered to be small |13j . In- 
deed, we have corroborated that the inclusion of an outer 
isothermal sphere does not affect the conclusions of this 
paper. 

With the use of the expression for the projected mass, 
A%fdm(£*)j the lens equation simplifies to 
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FIG. 1. The lens equation of a SFDM halo model, Eq. (jSJ, as 
a function of A. The roots define the Einstein radius, #*e, and 
its local maximum (minimum) the critical impact parameter, 
/3* C r- Both quantities are well defined only for values of A > 
A cr ~ 0.27, which is the threshold value for strong lensing. 



where to(£*) = M S fdm (£* ) /Pc^max i s a normalized mass 
function, evaluated numerically. Here /3* = DohP/r max 
and 0* = I?oL0/7'max are the normalized angular posi- 
tions of the source and images, respectively, and the pa- 
rameter A is given by 
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In order to avoid confusion with the self-interaction term, 
A, we have introduced a bar in the new parameter A. We 
have also defined the reduced angular distances dA = 
DaHq/c, and considered H = 100/i(km/s)/Mpc as the 
Hubble constant today, with h = 0.710 ± 0.025 [14] . 

In Fig. [l] we show the behavior of the lens equation |5]) 
as the A parameter varies (i.e. for different values of the 
combination p c r max ). Some notes are in turn: i) Strong 
lensing can be produced only for configurations with A > 
A cr ~ 0.27, and ii) For these configurations, only those 
with an impact parameter 1/3*1 < /3* CI can produce three 
images (note that the actual value of /3* cr depends on the 
parameter A, /3* cr (A)). 

These conditions on the SFDM profile are very similar 
to those obtained for the Burkert model in [15] ; this is not 
surprising because both of them have a core in radius. In 
that sense SFDM halos are analogous to those proposed 
by Burkert |16j . but with the advantage that their prop- 
erties are clearly connected to physical parameters in the 
model. 

In Fig. [2] we show the magnitude of the Einstein ra- 
dius, 0*e, as a function of the parameter A, where for 
comparison we have also plotted the same quantity for 
the NFW [T7] and Burkert [T5] profiles. The minimum 
value of A needed to produce multiple images is higher for 
a SFDM halo, A£ FW = < A^ r urkcrt = 2/ir 2 < ~X^ DM ~ 
0.27. (Notice that there is an extra factor of 1/4tt in 
our definition of A when compared to that reported in 
Ref. [T5].) SFDM halos seem to require larger values of 
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FIG. 2. The Einstein radius, #*e, as function of A', for 
SFDM (solid line), NFW (dashed line), and Burkert (dotted 
line) halo models. Einstein rings of similar magnitude require 

pFW ^ -^Burkert ^ ^SFDM 



A in order to produce Einstein rings of similar magnitude 
to those obtained for the other profiles, but this is in part 
due to projection effects that have not been considered 
in this paper [13l [18] . 



III. LENSING VS DYNAMICS 

Taking into account that in SFDM models there is a 
critical value for the parameter A, A cr ~ 0.27, and con- 
sidering the definition in Eq. (5b I, we can write the con- 



dition to produce strong lensing in the form 
<[M 0P c- 2 ] > 473.68/!^ , 



PcTn 



(6) 



with /di S t = dos/dohdhs a distance factor. 

In order to evaluate the right-hand-side (r.h.s.) of 
Eq. we consider two surveys of multiply-imaged sys- 
tems, the CASTLES [TH] and the SLACS [20]. From 
them we select only those elements for which the red- 
shifts of the source and the lens have been determined 
(which amounts to approximately 60 elements in each 
survey), and calculate the corresponding distance fac- 
tor /dist for every element in the reduced sample. In 
CASTLES (SLACS) the distance factors are in the inter- 
val 4 < /di S t < 27, (6 < /dist ^5 25), with a mean value 
of /dist - 7, (/ dist ^ 11), and then the r.h.s. of Eq. M 
takes on values in the range 1400 - 9000, (2000 - 8500). 
Some representative elements from SLACS are shown in 
Table [I] (galaxy lensing). In terms of the mean values, 
the inequality in Eq. (pi) translates into 

p c r max [M pc- 2 ] > 2000 , (CASTLES) (7a) 
p c r max [M pc- 2 ]>4OOO, (SLACS) (7b) 

These numbers are an order of magnitude greater 
than those obtained from dwarf galaxies dynamics, 
Pc?"max[-M©pc~ 2 ] — 100, when interpreted using the same 
density profile [7]; see again Table [I] This is the main re- 
sult of the paper. Remember that the value of r max is re- 
lated to the fundamental parameters of the model, which 



are the mass of the scalar particle and the self-interaction 
term, and it remains constant throughout the formation 
of cosmic structure. 

We must recall that inequalities in Eq. ^ do not take 
into account the presence of baryons in galaxies. Gravity 
does not distinguish between luminous and dark matter; 
then the contribution of the former to the lens could be 
significant in some cases. For instance, for those systems 
in SLACS the stellar mass fraction within the Einstein 
radius is 0.4, on average, with a scatter of 0.1 [21"] . 

We have corroborated that our estimates in Eq. (|7|) are 
not sensitive to the inclusion of a baryonic component. 
To see that we add the contribution of a de Vacouler 
surface brightness profile [22] to the lens equation, 

M0,)=0,-X^-X lm J-^l. (8a) 

Here Ai um is a parameter analogous to that given in 
Eqd5bk, 
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and f(x) a dimensionless projected stellar mass, 
1 



(8b) 
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-2520g 2 + 5040g - 5040) + 5040] , (8c) 

with q = —7.76x~ 1 / 4 . The mass-to-light ratio, M/L 
(from a Chabrier initial mass function), and the effec- 
tive radius, r e , for each system in SLACS are reported 
in Ref . [2T] . With the use of Eq. Q strong lensing is al- 
ways possible. Consequently, we must impose a different 
condition to constrain the product p c r max in each galaxy, 
such as demand the formation of Einstein rings of certain 
radius. 

To proceed we use a small subsample of SLACS that in- 
cludes configurations with the minimum, maximum, and 
mean Einstein radius, and stellar surface mass density, 
respectively. This is because the new lens equation is a 
function of the ratio r eif — r e /r max ] then, to compute the 
magnitudes of the Einstein radii, we shall fix the value of 
r max a priori. 

Using the new lens equation we find the value of A that 
produces the appropriate Einstein radius for each of the 
elements in the subsample. This is done using two differ- 
ent values of r max : 5 and lOKpc. The resultant products 
PcTmax are compatible (in order of magnitude) with the 
inequalities obtained from Eq. ([6|. Only for those sys- 
tems in the subsample with a high stellar surface mass 
density the value of p c r max can decrease substantially, 
but it is important to have in mind that all the possible 
uncertainties associated to the distribution of the lumi- 
nous matter, like the choice of the stellar initial mass 
function, will be more relevant in such cases. In gen- 
eral, these estimations are sensitive to the details of the 
particular configuration, and a more exhaustive analy- 
sis, considering the complete sample, will be presented 
elsewhere. 
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DYNAMICS OF GALAXIES 


GALAXY LENSING 


Galaxy 


p c r m ax[M pc"' i 


Galaxy 


/dist 


PcTmax[M pC~' ! 


Ho II 


36.19 


J0008-0004 


6.61 


2029.68 


DDO 154 


66.47 


J1250+0523 


8.46 


2832.41 


DDO 53 


67.53 


J2341+0000 


9.12 


3053.38 


IC2574 


81.89 


J1538+5817 


11.74 


3930.44 


NGC2366 


85.45 


J0216-0813 


13.03 


4362.44 


Ursa Minor 


104.72 


J1106+5228 


15.74 


5269.75 


Ho I 


120.23 


J2321-0939 


16.23 


5433.80 


DDO 39 


145.94 


J1420+6019 


19.72 


6602.26 


M81 dwB 


265.87 


J0044+0113 


25.26 


8457.05 



TABLE I. Estimates of the product p c r max for different galaxies. Left. As reported in Refs. [7], using galactic dynamics. 
Right. Derived from equation Q in this paper; recall that these values represent a lower limit (here we show only a representative 
subsample of the SLACS survey). Note the difference of an order of magnitude between the values of p c r ma x for dwarf galaxies 
in the local universe, and the lower limit of this same quantity for galaxies producing strong lensing at z ~ 0.5. 



IV. DISCUSSION AND FINAL REMARKS 

We have shown that a discrepancy between lensing and 
dynamical studies appears if we consider that the SFDM 
mass density profile in Eq. ([!]) describes the inner regions 
of galactic halos at different rcdshifts, up to radii of order 
5 — lOKpc. More specifically, we have found that lens 
galaxies at z ~ 0.5, if correctly described by a SFDM halo 
profile, should be denser than dwarf spheroidals in the 
local universe, in order to satisfy the conditions necessary 
to produce strong lensing. 

In principle nothing guarantees that halos of differ- 
ent kind of galaxies share the same physical properties. 
Our studies took into account galaxies that are intrin- 
sically different in terms of their total mass and baryon 
concentration. While dwarf galaxies show low stellar sur- 
face brightness, stellar component in massive, early type 
galaxies is typically compact and dense. 

In the standard cosmological model the evolution of 
DM halos may trigger differences in concentrations for 
halos with different masses due to differences in the as- 
sembling epoch; smaller halos collapsed in an earlier and 
denser universe, therefore they are expected to be more 
concentrated. However, it is also well known that the 
presence of baryons during the assembly of galaxies can 
alter the density profile of the host halos and modify 
this tendency, making them shallower (supernova feed- 
back [23])) or even cuspier (adiabatic contraction 24|). 
Therefore, the stellar distribution may reveal different 
dynamical evolution for low and high mass halos trig- 
gered by galaxy formation. 

For SFDM, the dynamical interaction between baryons 
and the scalar field may also modify the internal halo 
structure predicted by the model, Eq. ([I]), clarifying the 
discrepancy. For instance, if the concentration of stel- 
lar distribution were correlated with that of the halo, 
like in the adiabatic contraction model when applied to 
standard DM halos, this may explain our findings. But 



at this time it is unknown how compressible SFDM ha- 
los are, and if such effect will be enough to explain our 
results, because there are no predictions on its magni- 
tude. If the modification triggered by baryons were insuf- 
ficient, then it might be suggesting an intrinsic evolution 
of SFDM halos across cosmic time. For example, if big 
galaxies emerge as the result of the collision of smaller 
ones, then the central densities of the resultant galaxies 
would be naturally higher; after all, r max is a constant 
in the model, and one would expect that total mass is 
preserved in galaxy-galaxy mergers. At this point we do 
not know which of these two mechanisms, the intrinsic to 
the model, or that due to the evolution of SFDM halos 
in the presence of baryons, is the dominant one. In that 
sense, a theoretical description of these processes may be 
very useful and welcome. 

A full picture requires a distribution of values for the 
central density generated from the evolution of the spec- 
trum of primordial density perturbations after inflation. 
Such a result is not available now, but it is possible to 
start tracing this distribution with galaxy observations. 
We present, for the first time, observational constraints 
on the dynamical evolution of SFDM halos in the pres- 
ence of baryons, that must be considered for future semi- 
analytical/numerical studies of galaxy formation. 
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